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I, INTRODUCTION 


A well-known programming problem is the traveling-salesman 
problem. As classically stated, a salesman must start from his home- 
town, travel to each of n-l other cities, and return home. He is re- 
quired to enter each of the n cities exactly once and leave each city 
exactly once, with the understanding that his entire trip is one 
continuous circuit. There is no requirement that the distances be- 
tween cities be symmetric, that is, the distance from city i to city j 
need not be the same as that from j to i. The salesman's objective is 
to select a sequence of cities in such a way as to minimize the total 
distance traveled, 

There are several algorithms which solve the traveling-salesman 
problem, but none offers the last word in terms of speed and ease of 
computation. The purpose of this investigation was to examine a newly 
proposed solution procedure developed by Professor Harold Ca 
His algorithm looked very promising in terms of hand computation on 
small problems, but this study focused on its desirability for 
handling larger problems, such as the forty-city case. The greatest 


effort went into programming the algorithm for the computer, then 


letting the computer serve as the primary tool of analysis. 


lcreenberg, Harold, Professor, Department of Operations Analysis, 
Naval Postgraduate School, Monterey, California 


A, DEFINITIONS 
The distance (cost) between city i and city j is denoted Ci; for 


i = L,...,neand j = Li ..4m. lagcencra ace 


ij need not equal eae The 


decision variable se equals one if the salesman is to travel directly 
from city i to city j and equals zero, otherwise. Since traveling 
directly from city i to city i is not meaningful, Xs must always 
equal zero for all i. Most solution procedures meet this requirement 
by assigning arbitrarily large values to Cae feo qamerc el A GH 

Associated with each city is the number of steps required by a 
solution to arrive at that city. For example, if the salesman is to 
go from city A to B to C and back to A, he arrives at B in one step, C 
in two steps, and back to A in three steps, Unless otherwise stated, 
the following discussion takes city one to be the city of origin and 
city n to be the last city visited prior to returning to city one, 

A solution must provide a tour in order to be a traveling-salesman 
solution. A tour is defined as a sequence of transitions from city 
one through all other cities and back to city one, such that each city 
is entered exactly once and left exactly once, and it must take exactly 
n steps to return to city one. The latter clause guarantees there are 
no discontinuities in a tour, Deleting this clause yields the defini- 
tion of a subtour. It is possible for a candidate solution to contain 
as many as n/2 subtours and still meet the requirement that each city 
is entered and left exactly once. To illustrate the concept of tours 
and subtours, consider these solutions to the four city problem: city 
1 to 3 to 4 to 2 to 1 provides a tour; whereas city | to 2 to l@aue 


city 3 to 4 to 3 provides two subtours. 


B. ALTERNATE FORMULATIONS 

The traveling-salesman problem is essentially a combinatorial 
problem and there is no one way to solve it. The most direct approach 
is to select a starting point, evaluate the distance (cost) over each 
of the (n-1)! possible tours, and then choose the tour(s) with minimum 
total cost. A dynamic-programmjing formulation reduces the number of 
tours one must consider but the tables required can quickly exceed 
computer storage capacity. Branch-and-bound techniques can also be 
employed but suffer from similar limitations as dynamic programming. 
A fourth possibility might be to construct a solution procedure using 
an algorithm for finding minimal-cost flows in a network. 


The traveling-salesman problem can be stated as an integer program- 


ming problem. Hadley / Ref. 6/7 and Dantzig /Ref. 1/ present essential- 


ly the same formulation. Dantzig lets Rik = 1 or O according to 


t 


whether or not the salesman travels from city i to city j on the t 


he 


3 


Seep, where ijyj,itt =a7,...jn. Mien, defining tas | = ea 


states the problems as follows: 


minimize = oy Se = Z 
L,J,¢ a tae 


subject ‘to a = oe G6 tails ee 
bhai k 4,k,t+1 

and <a” = Ger —*1 ,..... 
Ie 1 see 

and ne = or lefer abl i,i.t. 


The first set of constraints requires that if the salesman arrives at 
city j in t steps, he must leave that city on step t+l. The second 
sets requires that the salesman leave each city i once and only once. 
: 3 , 2 ' , 
Solving this system of n + 1 variables and n + n equations with an 


integer programming algorithm yields a minimum-distance tour. 


The integer programming approach has several shortcomings. The 
third subscript t on x greatly increases the number of variables, and 
the number of constraints in the formulation is very large. Further- 
more, integer programming algorithms normally require many more itera- 
tions than ordinary linear programming techniques. 

There is a general class of indirect solution procedures that 
takes advantage of a special feature of the traveling-salesman problem. 
If one temporarily ignores the fact that a solution must provide a 
tour, the salesman's task can be formulated as a simple assignment 
problem, perhaps the most readily soluble of all linear programs. 

The assignment formulation is achieved by letting i = 1 orad 
according to whether or not the salesman is to travel from city i to 


je Let Ges be arbitarily large for all i. Then, 


i ty 


minimize 2 = c x 
1 ) 


pea 


subject to i Ki == oh Cj). trleyeeeegn) 
and aes x =l Gin = do Sa 
J 1) 


Notice that there are only aa variables and 2n constraints. One 
constraint may be dropped because exactly 2n-1 are linearly independent. 
Furthermore, since the assignment problem is a form of the transporta- 
tion problem, and since all ee ae expressions equal one, any 
solution should yield 0 or 1 values for the xi, 8 dmePet . elif. Thus, 
there is no need to add the integer constraint to the formulation 


above. Possible solution procedures include the simplex method, the 


familiar transportation tableau, or an algorithm for finding minimal- 
cost flow in a network / Ref. 4/. 

An assignment outcome does not necessarily provide a traveling- 
salesman solution, but it does satisfy the requirements that each 
city is entered exactly once and departed exactly once and that a. = 0 
for all i. Clearly, since the salesman's problem includes the added 
constraint that any set of values for the Xi 4/8 must provide a tour, 
his minimum-distance solution must be at least as great as the minimum 
assignment distance. 

The so-called indirect procedures, alluded to above, typically use 
an assignment solution as a starting point. The optimal value of the 
assignment objective function provides a lower bound to any traveling- 
salesman solution. A non-optimal assignment solution may still provide 
a lower bound while being closer in magnitude to the salesman's minimum 
total distance. In the sense that it is closer, the latter lower 
bound is said to be "better" than the former. In practice, however, 
it is difficult to find a value larger than the optimal assignment 
solution and prove that it is, in fact, a lower bound to the traveling 
salesman problem. 

Professor Greenberg's algorithm uses the optimal assignment solu- 
tion as its point of departure because it provides a guaranteed lower 
bound and because special use is made of the coefficient vectors in 
the final simplex tableau. After the solution of the initial assign- 
ment problem, the algorithm simply manipulates the coefficient 
vectors until a minimum-cost tour is achieved. This procedure is 


explained in detail in the next chapter. 


W. L. Eastman also begins with the optimal assignment solution and 
if it should provide a tour, then the problem is solved. If, in the 
more likely event, it does not, he selects the subtour with the least 
number of variables and sets each equal to zero, one at a time. For 
each areas he resolves the original assignment problem but with the 
added constraint. In practice this can be achieved by setting 
Cis = ©© corresponding to the Xi required to be zero. If a tour is 
not obtained, he branches from that assignment problem with the 
smallest-valued objective function. Again the least-membered subtour 
is chosen to create new assignment problems, and so forth, until a 
tour is achieved. /Ref. 8/ 

J. D. C. Little and others have devised a branch-and-bound 
algorithm that starts with a lower bound, not to the traveling salesman 
problem, but to the assignment problem. Instead of obtaining the 
optimal solution to the assignment problem at each step, the smallest 
element is subtracted from each row of the cost matrix, and then the 
smallest element from each column of the result. The sum of the re- 
maining elements is said to be a valid lower bound on the optimal 
assignment solution, At each branching step two new problems are 
created corresponding to xij = O and ie = 1. The ae chosen for each 
branching step is the one that yields as large a bound as possible 
when set equal to zero. Branching proceeds until a tour is obtained. 
Little's method normally requires many more branches than does Eastman's, 
but there is less computation at each step, Some performance results 
are available on Little's algorithm. The mean execution time on an 
IBM 7090 computer for 100 thirty-city problems was 58.5 seconds and the 


™ 


average for five forty-city problems was 8.3/7 minutes. /Ref. 9/ 
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C. PLAN OF THE REPORT 

Chapter II presents a detailed description of the solution procedure 
suggested by Professor Greenberg. His proposed algorithm is complete 
in the sense that the criteria for proceeding from one step to the 
next are completely specified and guarantee an optimal solution. In 
some instances, however, it is not readily apparent how one can most 
efficiently meet these criteria. Chapter III deals with these ques- 
tions of technique as it describes how the algorithm was programmed 
for the computer. When there were alternate computational approaches 
possible, an attempt was made to find a good one. There is no claim 
that the particular techniques devised to meet the algorithm criteria 
are the ones Professor Greenberg, himself, would have chosen. 
Chapter IV outlines the salient results of computational experience 
with the computer program, The algorithm and the program are evaluated 


in the light of this experience. 


bl 


LI. THE ALGORITHM 


An algorithm is defined by Webster as "a rule or procedure for 
solving a mathematical problem that frequently involves repetition of 
an operation.'' This chapter seeks to present Professor Greenberg's 
algorithm for solving the traveling-salesman problem in the logical 
sequence of its development, concluding with a concise statement of 


the algorithm, itself. 


A. THE TRAVELING-SALESMAN PROBLEM 

A traveling salesman is faced with the problem of finding the 
shortest route (or the least-cost route) through n cities. Further- 
more, he must meet the following conditions: 1) he must enter each 
city once and only once, 2) he must depart from each city once and 
only once, and 3) his route must be a tour, that is, there must be no 
discontinuities in his route. There are (n-1)! possible solutions to 
this problem, since the salesman can go to any of n-1 cities on his 
first step, then to any of the remaining n-2 cities on his second 
step, and so forth, until he returns to his starting point on the ine 
step. Optimal solutions, of course, are found within this set of 
(n-1)! tours. 

Clearly, the above is merely an allegorical statement of a problem 
that has relevance to a variety of situations. For example, suppose 
fixed electronic components must be connected in a continuous circuit 
such that the least amount of wire is required. Disregard of condi- 
tions (1) and (2) might leave some components out of the circuit 


entirely or, on the other hand, cause short circuits. Violation of 


le 


condition (3) would result in two or more subcircuits, with current 


flowing in only one. 


B. THE ASSIGNMENT FORMULATION 

Modifying the traveling-salesman problem by disregarding the third 
requirement above, yields a problem which is readily adaptable to the 
assignment formulation, perhaps the most easily soluble of all program- 
ming problems. The assignment problem is usually stated in terms of 
assigning n men to n jobs in such a way that overall productivity is 
maximized, assuming that a productivity index can be associated with 
each man-job combination. There are n! possible solutions to this 
problem, since the first man can be assigned to any of the n jobs, the 
second man to any of the remaining n-1l jobs, and so forth. 

Dropping the requirement, then, that a solution must provide a 
tour, the modified salesman problem is to assign n cities to n cities 
(the same n cities) such that the sum of city-to-city distances is 
minimized. The possibility of assigning city i to itself is precluded 
by giving each C;, an arbitrarily large value. The possibility of 
other types of subtours appearing in a solution, however, cannot be 
restricted without destroying the assignment formulation. 

Since no city is ever assigned to itself in the modified salesman 
problem, the size of the solution set is smaller than n!, the usual 
size for assignment problems. The table below illustrates solution- 


set sizes for several selected problems: 


rS 


Number Number of Number of Number of 


of cities Assignment Modified Traveling- Solutions that 
(n) Solutions Salesman Solutions are Tours 
(n!) ne) = noe (n-1)! 
2 2 i, i 
3 6 2 Z 
4 24 2 6 
5 120 44 24 


A generalized expression for the third column above was not derived, 
but it is clear that the fraction of the time a modified salesman 
solution is also a tour is greater than (n-1)!/n! = 1/n. 

An optimal solution to the modified salesman problem provides a 
lower bound on the minimum-distance solution to the original problem. 
This is true, because the modified problem differs from the original 
only in that it has fewer constraints, If the optimal assignment 
solution to the modified problem provides a tour, then the original 
problem is solved; if not, further steps must be taken, 


To express the assignment formulation in the notation of Chapter I, 


let a 1 or O according to whether or not the salesman is to travel 
from soit) Woy Cig | itt 8 ellie Slee case, ieee o be the distance between 
i} 
CACYselueand jeebotea | led ponte, GEE Cie arbitrarily, dargesfior alias 
Then, 
minimize = =. ate gtk 
. ae 
subject to aa et = ery 
i ij 
and ae See | ile, ee 
i i j 
and x, BeO-tor all 14 
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Alternately, this formulation can be expressed in matrix notation 


as follows: 


Z 
minimize Cum = 2 
subject to Ax = b 
and x20 


where the definitions of co! x,b, and A can be indicated most con- 


veniently by example. Consider the three-city problem: 


iN 
C= (Cy1, SyQ2 C439 Cny> Cone £239 $312 ©3929 %33)> 


X11 1 
X19 i 
X13 1 
Xo L 
Xx - Xoo ‘ b = 1 ‘ 
X53 1 
X39] k 
X39 I 
X33 1 
1 I 0 0 0 0 0 
0 0 1 1 1 0 0 
0 0 0 0 0 1 z 
and A = 0 0 1 0 0 1 0 
1 0 0 1 0 0 1 
0 1 0 0 i 0 0 





i ae ..42 a. 2 
In generals e" ¥o 1 x n*,”x is n® x l,°bD Ys mn” % FF, and A is 2m x me. 


th X,. Gi m 


Notice that since the ye column of A corresponds to the k ij 
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the aa element and the (aii element of the Pa column of A equal 
one and all other elements in that column equal zero, 

The assignment problem has several important features that are 
listed below without proof. (see references 3 and 5 for proofs.) 

1. The assignment problem has a feasible solution. 

2. The rank of the matrix A is 2n-l, thus any row may be deleted 
as redundant. 

3. Any set of 2n-1 linearly independent column vectors from A 
forms a basis B. 

4, The inverse of any basis B is composed entirely of zero, one 
and minus-one elements. 

5. Any basic feasible solution to Ax = b is degenerate. In 
fact, there are exactly n basic variables equal to one and exactly 
n-l equal to zero, assuming all nonbasic variables are arbi tran die 


set equal to zero. 


C. THE SIMPLEX SOLUTION OF THE ASSIGNMENT FORMULATION 

There are alternate methods for solving the assignment problem but 
the following discussion assumes some variant of the simplex method is 
chosen. The final simplex tableau contains not only the optimal solu- 
tion to the assignment problem, but other interesting information, as 
well. Before proceeding, a few more notational conventions and 
definitions are stated below: 

1. Redefine the coefficient matrix A sicr, that one row (any row) 
is deleted. 

2. Temporarily redefine the vector x such that the subscripts on 


; z 
its elements run from 1 to n° rather than 1,1 to n,n. 
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3. Make the same change in subscript notation for the elements 
of ~ 


.th 2 
4. Denote the j column of A as af ner all j= lye 6 


5. Let B denote a basis in A, that is, any set of 2n-1 column 


vectorsS,;yin A. 


6. “Let Xp denote the vector of XS corresponding to the a,'s in 
B. Call the elements of Xp the values of the basic variables. Note, 
since Xp = like Xp is simply the vector obtained by summing each row 
of B. Requiring Xpz 0, then, insures its elements are ones and zeros, 


exclusively. 


T 
ve Het Ce denote the vector of c.'s corresponding to the elements 


in Xp. 
-1 
8. Define ys = B ae Note that the y,'s are composed entirely 
J 


of zero, one, and minus-one elements. 
9. Define the relative cost for x. as ¢, = ce -cly, 
j j , 

The simplex procedure involves finding a basis B such that x_ = 
Fw minimizes z. A feasible solution is recognized when the elements 
of xp are zeros and ones, exclusively, and an optimal solution, when 
~~ = 2 
oy a 0 for all J = L 2a Sen e 

Notice that c, = 0 always for all ae in Xp because for such 


T 7 
er ee Cp y;. If — = 0 at optimality for a nonbasic Xi then an 


J J J 
alternate optimum can be obtained by bringing that Xs into the basis. 
In this case, the simplex methodology specifies which basic variable 
is to be replaced in order to maintain feasibility. 


To appreciate the impact of the relative costs on the value of z, 


consider the following relationship: 


J 


Q=e + = C.X 
jer 3 
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where és is the new value of zg and R is the index set of all nonbasic 
variables. This equation shows that if nonbasic variables are set 
equal to one, z is increased by the sum of their relative costs. 

Of course, if one or more nonbasic variables are set equal to one, 
some change must take place in the xp vector. This is true since an 
assignment solution requires exactly n variables equal to one, includ- 
ing the nonbasic variables. In order to measure such effects on Xp-s 


consider the following equality: 


m= %,- See YE, 


B B j€R jj 


A 
where’ x, is the new set of values for Xp Normally a = QO for all j in 


R, but if one or more such nonbasic variables are set equal to one, the 
relation above changes the values of the basic variables. Setting one 
or more x's not in Xp equal to one is said to be constraint feasible 


A , 
as long as the elements of xX, are zeros and ones, exclusively. 


D. STATEMENT OF THE ALGORITHM 

The algorithm under study begins with the optimal assignment 
50 GENO to the modified traveling-salesman problem, If the optimal 
assignment solution provides a tour, then the original problem is 
solved, Otherwise, the algorithm begins to introduce combinations of 
nonbasic variables at the one level until a tour is achieved. 
Complexities arise in the measures taken to assure that 
combinations are considered in ascending order of their relative costs. 
If these steps are not taken, there is a possibility that nonbasic 
variables not yet tried yield lower-cost tours. 

In the light of the background developed in the previous sections, 


the algorithm is as follows: 
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1. Drop the constraint that a traveling-salesman solution must 
provide a tour. 

2. Solve this modified problem as an assignment problem using a 
variant of the simplex method. 

3. Check the optimal assignment solution for a tour. If a tour 
is provided, stop. Otherwise, continue. 

4. Drop simplex tableau vectors corresponding to the original Xs 5 
variables. Drop tableau vectors corresponding to basic variables. 
There are at -3nt1 vectors remaining. Let m = mee3nti. 

5. Rank tableau vectors in ascending order of their respective 
relative costs. 

6. Let xi = =; = gs and Se = c. where i = 1 if c, is first 
in the ranking, i = 2 if c. is second, and so forth, through i =m if 


J 

C; is last in the ranking. (Retain information required to reconvert 
subscripts as needed when checking for tours.) 

7. Let i=0. 

8. Increment i by one. 

~~ TS x; = 1 is not constraint feasible and i = 1, return to 
geep , Tf x. = | Ts wet “constraamt feasible and 171, goto sicep 1. 
If x; = 1 is constraint feasible, continue. 

20. WE x. = 1 and c ae ae provide a tour, go to step 21. Other- 
wise, return to step 8 if i = 1; continue if i >? 1. 

11. Generate all combinations of x, = 1 With = = l, a Se. 
X,_, = 1. (Note this is construed to mean x, = 1 appears in each 


combination along with one or more of the other variables indicated set 


i-1 ; 
equal to one. There are 2 - 1 such combinations.) 
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12. Disregard all combinations with relative costs greater than 


QO 
Le) >) 


(Define Ctl = 
13. Disregard all remaining combinations that are not constraint 
feasible, 
14, Check remaining combinations, in ascending order of their 
relative costs, for tours. As soon as a combination is found to pro- 


vide a tour, go to step 21. If no combination yields a tour, continue. 


Losec £. isey2. and/ orwes 


i+] = Ci» return to step 8. Otherwise, 


continue. 
LO, ubet jkeaas 
17. Decrement k by one. 


18. Generate all combinations of x, = Lawith x =~ IL: Xy = lowe 33 


1 


“ear © 

19. Same as steps 13, 14, and 15. 

20, If k > 2, returmeto step L7cpelt kise25 mebueneto Racers 

21. If coming from step 19, let i= k. Otherwise, do not change i. 

22. Compute total distance for the tour found. Letting K be the 
index set of nonbasic variables included in the solution, the total 
distance for this tour is £ = Zz 1, Ch. 

23. If i 2, go to step 32.4 Othemiise Mconernuec, 

24. If all combinations of nonbasic variables preceding Xs have 
already been checked against the current value of 2-7 as an upper 
bound, go to step 32. Otherwise, continue. 


25. Decrement i by one. 


26. Generate all combinations of x = 1 with x, = l, X. ssl. ae 
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28. Disregard all remaining combinations that are not constraint 
feasible. 

29, Check remaining combinations, in ascending order of their 
relative costs, for tours. As soon as a combination is found to pro- 
vide a tour; go to step 30. If no combination yields a tour, go to 
step 31. ‘ 

30. Update the value of 9 for the latest tour found and continue; 

31. If i? 2, returti to step 25. Otherwise, continue. 

32. The latest tour found and its associated % is an optimal 
traveling-salesman solution. Stop. 

Since eventually all combinations of nonbasic variables are con- 
sidered, the algorithm must eventually find a tour. Suppose a tour 
is first found using x, = 1 in some combination with the preceding 
nonbasic variables. This tour's relative cost is less than or equal 
to Cay because only combinations meeting this criterion were considered, 
Thus, there is no need to try combinations beyond Xs There is a need, 
however, to reconsider combinations prior to Xa For example, x, = 1 
in combination with xX, = 1 may have been rejected earlier because its 


| 2 A 
relative cost is greater than c but it may be less than z-z, Thus, 


og 
proceeding backward through the nonbasic variables using the current 


value of Q=2 as an upper bound guarantees the minimum 4 will be found. 
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III. THE COMPUTER PROGRAM 


The algorithm was programmed for the IBM System/360 Model 67 


Operation System using FORTRAN IV (H level) programming language 


ef, ae A listing of the complete program is provided in Appendix 


l~ 


a 


Several routines used in the program were not developed as part 
of this study but were modified for use here as needed. Subroutines 
ASSIGN, LIP, ELEM, and PIVOT were supplied by Professor Greenberg and 
are used to solve the initial assignment problem. Subroutine RAND is 
a simple routine for generating city-to-city distances, uniformly 
distributed from 0 to 100. The distance from each city to itself is 
arbitrarily set equal to 999, Subroutine TIMEIT is used to measure 
algorithm execution times in seconds, Part of TIMEIT appears in 
assembly Ten cae 

The following list gives a brief description of the various 
routines developed especially for this study. 

1. The Main Program follows the sequence of algorithm steps as 
stated in Chapter II. Provision is made to either read in a cost 
matrix or to generate one randomly. Other provisions are made to 
abort if something unexpected goes wrong and to begin execution on 
the next problem. 

2. Subroutine PRINTC writes out the cost matrix. Also, alternate 
storage space is provided for the costs because they get overlayed 
in the assignment routines by the relative costs. 

3. Subroutine PRINTA writes the optimal assignment solution in 


tableau form. If the problem is larger than seven cities they ys; 
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column vectors are not printed because there are too many to fit ona 
page. Subroutine PRINTA is also used again to print the optimal assign- 
ment solution after artificial variables remaining in the basis have 
been replaced, 

4. Subroutine SOLCHK checks constraint feasible solutions for 
tours. If a tour is found, it.is printed. 

5. Subroutine RANKC ranks the relative costs in ascending order of 
magnitude. Actually, the costs are not moved around in storage, but 
second-order subscripts are assigned to indicate the desired ranking. 

6. Subroutine REPART was devised to replace any artificial 
variables left in the assignment solution and to select a redundant 
row for deletion from the final tableau. 

7. Subroutine LINCOM is the most important subroutine in terms 
of the algorithm at hand. This subroutine generates the combinations 
of the nonbasic variables, compares the relative cost of each combi- 


nation with an upper bound, and checks for constraint feasibility. 


A, DEFINITIONS 

Three common blocks (storage areas) are used in the program. 
Common Block ASGBLK is used exclusively with the four assignment 
routines. Common Block SHARE is used to pass information from the 
assignment routines and from Subroutine RAND to the traveling-sales- 
man algorithm routines. Common Block TSBLK is used to pass values 
among the algorithm routines. For easy reference, this section 
defines the variables in common in the order of their appearance. 
Variables in ASGBLK are not listed because they are only needed 
within the four assignment routines which are not part of this 


study. 
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1. Common Block TSBLK 

a. COST is the sum of the relative costs included in a 
particular combination of nonbasic variables, each taken at the one 
level. 

b. CSTLC(1600) is initially used as alternate storage for 
the city-to-city costs. After artificial variables have been replaced, 
it is used to store relative costs of combinations in ascending order 
of magnitude, 

c. NLC(1000) is used to store a combination identifying 
number corresponding to each value in CSTLC(1600). With these numbers, 
combinations can be regenerated as needed, 

d. 1C(1600) is used to store the subscripts of the variables 
associated with the ranked, relative costs. For example, the first 
value in IC(1600) is the subscript of the nonbasic variable with the 
smallest relative cost. 

e. KC(40) is used in Subroutine PRINTIC to write out the city- 
to-city costs as integers, 

f. KB(120) is used to store the new values of the basic 
variables when one or more nonbasic variables are set equal to one. 

The one-values for the nonbasic variables are augmented to this vector 
starting after the last basic-variable value, (Also, see description 
of IB(120) below.) 

g. IAX(40) is used in Subroutine PRINTA to write out a row 
of the final assignment tableau. 

h. IROW(40) and ICOL(40) are used in Subroutine SOLCHK to 
store the first and second subscripts, respectively, of variables after 


these have been reconverted from single-subscript form. 
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i. LC(40) is used to store a binary number up to 40 digits in 


length. If the eh digit is d-one, then the yeh ranked, nonbasic 


variable is to be included in the particular combination at the one 


o ranked, nonbasic variable 


is not to be included, LC(40) can be used to represent up to gee - l 


th 
level. If the k digit is zero, then the k 


different combinations, one at a time. 

j. KZ is the truncated value of Z in Subroutine PRINTA. In 
Subroutine SOLCHK, KZ is the truncated vailte OL" XZ =" Z + eter. “Fr the 
city-to-city costs are integers, KZ = XZ always. 

k. NC = N-3MM+1 is the number of nonbasic variables, excluding 
Kes variables and artificials, 

1. J is used in the Main Program to index the progression 


through the ranked, nonbasic variables. 


m. @¢3 = IC(J). 


a Lie DECC, (CS) wand 82 = Bee (2, ICH). 

o. Jl = J-1. 

p. CST is the current upper bound against which combinations 
are compared, Any combination of nonbasic variables at the one level 
with a relative cost greater than CST is disregarded. 

q. XZ = Z4+COST. Subroutine SOLCHK truncates this value when 
writing it out. (See definition (j) above.) 

2. Common Block SHARE 

a. C€(1600) is used initially to store the city-to-city costs. 

After the assignment solution, the relative costs are stored here. 


(Also, see description of CSTLC (1600) above.) 


b. B(80,80) is the inverse of the optimal assignment basis. 
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c. H(80) is used to store the values of the basic variables 
in the optimal assignment solution. (Also, see description of KB(120) 
above.) 

d. 2Z is the optimal assignment total cost. 

e. LEG(2,1600) is a matrix of values used to convert single- 
subscripted variables back to their original double-subscripted form. 


For example, x, = X.. if i = LEG(1,k) and j = LEG(2,k) - MM where MM 


1J 
is the number of cities. Furthermore, LEG provides a convenient device 
-1 
for obtaining the product B ay = Yue It is known that the je and 


th 
( j+MM) elements of a, are ones, if i and j are defined as above, and 


all other elements of a, are zeros. Thus, rather than multiplying 


= th 
B : times ays simply add the aos and (j+MM) elements of each row of 


a to produce the vector Ye 

f. 1IB(120) is used to store the subscripts of the optimal 
assignment basic variables. Negative subscripts indicate artificial 
variables. Starting just after the last basic-variable subscript, 
nonbasic-variable subscripts are appended to this vector for each 
nonbasic variable set equal to one in a particular combination. (Also, 
see description of KB(120) above.) 

eg. MM is the size of the problem, that is, the number of 
cities, 

h. M is the number of constraining equations in the assignment 
formulation. Initially, M = 2(MM), but in Subroutine REPART it is 
reduced by one, 


z 
i. N = (MM) is the total number of variables, excluding 


arta fa-cialiscs 
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B. DESCRIPTION OF SELECTED SUBROUTINES 
The Main —s was written to conform closely to the algorithm 
as presented in Chapter II and can be followed readily from the listing 
in Appendix A, Further amplification of Subroutines LINCOM, SOLCHK, 
and REPART is considered necessary, however, in order to clarify the 
workings of the overall program. Subroutine REPART was devised only 
because-the assignment routines used leave artificials in the basis 
and do not drop a redundant equation, Subroutines SOLCHK and LINCOM 
represent techniques created to fulfill algorithm requirements and, as 
such, are the most important products of this study. 
-1l. Subroutine REPART 

Subroutine REPART scans the subscripts of the basic variables 
returned from the assignment routines. As soon as a negative sub- 
script is an its corresponding artificial variable is replaced 
in the basis by_.the first, ranked, nonbasic variable that has a nonzero 
element in the ee row of its Y; vector. Using this nonzero 
element as a pivot element, standard row transformations are performed. 
The result is that the nonbasic variable is brought into the basis and 
the artificial is dropped. The optimal value of the objective function 
does not change because the artificial variable was necessarily at the 
zero level implying its replacement must also be at the zero level. 
After replacing an artificial, relative costs for the remaining non- 
basic variables are recomputed and these variables are then ranked in 
ascending order of their new relative costs. The above process is 
repeated until all artificials have been removed from the basis. 

It is known that one of the assignment-tableau rows is redundant. 


If, while trying to replace an artificial, a nonzero pivot element 


od 


cannot be found, that row is earmarked for deletion, Otherwise, after 
all artificials have been replaced, each row is scanned for all zero 
elements. If such a row is found, it is dropped. Otherwise, the last 
tableau row is arbitrarily selected for deletion. If other than the 
last row is dropped, remaining rows are shifted upward in the tableau 
eo fil leehe sap, 

Normally, control is returned to the statement following the 
calling statement in the Main Program. If for any reason, however, 
more than one row has gone to all zeros, control is returned to state- 
ment 10 in Main and execution on the next problem is started. 

2. Subroutine SOLCHK 

Subroutine SOLCHK examines an assignment solution to see if it 
provides a tour, Its approach is to count the number of steps required 
to return to the city of origin. If this number is less than the total 
number of cities MM, then subtours exists and control is returned to 
the statement following the calling statement in Main, If the number 
of steps required is exactly equal to MM, the solution is a tour. In 
the latter event, the solution is printed with its associated total 
distance and control is returned to the statement in Main indicated in 
the parenthesis of the calling statement, 

To illustrate the bookkeeping involved, consider the two fol- 


lowing four-city assignment solutions: 


steps 
Xg = Xo 4 = 1 T 
XQ = x7 3 = 1 

a 


28 


steps 


X19 = Xs 9 = 1 1 
X13 = X saa = 1 3 
% Xd = 1 Z 
Xy = X19 = 1 4 


Both are assignment solutions, but only the second provides a tour. 

Notice that the routine will keep scanning the KB vector until 
the amy &™ one-value is found and then scanning is terminated. This 
Bwiseck is justified because only constraint feasible solutions are 
sent to SOLCHK and these solutions necessarily have exactly MM 
variables equal to one, including any nonbasic variables set equal 
GO erre. 

3. Subroutine LINCOM 

The essence of the algorithm lies in altering the optimal 
assignment solution by setting selected nonbasic variables equal to 
one until a tour is achieved. Subroutine LINCOM was devised to generate 
all the required combinations of nonbasic variables. An attempt was 
made to program LINCOM so that the number of combinations necessary to 
consider be held to a minimum. The objective was to generate them in 
ascending order of total relative costs so that as soon as a combina- 
tion's cost exceeded the current upper bound, no further combinations 
need be considered, This objective was only partially achieved. 

The most difficult problem lay in the fact that there is no 
way to completely predetermine the ranking of all possible combinations. 


Consider the case where Sue COME. c= Ce and Xo = 1 is found not 
constraint feasible. Before testing a 1, all combinations of 


ty = 1 with the preceding nonbasic variables must be tried because some 


may have total relative costs less than or equal to Cee 


29 


One way to test the combinations would be to try x, = i al 
and xy, =.1 singly with Xo = 1, then in pairs, then triplets, and 
finally all four. A better approach, however, is to set up the four- 
digit binary number 0000 and then simply increment it by one, fifteen 
times. Number the digit positions from one through four, from right to 
left. Then, if a one appears in the i position where na Ll, sles 
let 2 = 1 appear in the combination. The possible combinations can 


now be represented as shown below: 


Identifying 
Number Combination 
1 0001 
Z 0010 
3 0011 
4 0100 
S: 0101 
6 0110 
7 Oe 
8 1000 
9 1001 
10 1010 
11 Low 
PZ 1100 
7, de IEIEONE 
14 1110 
15 1 lat 
Combination number 9, for example, represents Xe = 1 with x, = 1 and 


x, = 1 for a total relative cost of cs + oy, + Cy. 
Notice that the following generalizations can be made about 
the respective total relative costs of the combinations above: 
Fle #12 ee, elon 
#4 #5 = #6 = #7, 
#8 = #9 => #10 = #11, 


and tip <—®% 41398 eae: 


30 


Furthermore, #3 < #5 © #9, but #3's relationships with #4 and #8 
Cantiot be stated, in general, Also, #7 & #11 & #13, but #7's 
relationships with #8, #9, #10, and #12 are unknown. Finally, #11<4 
#13, but there is no way to tell in advance which is smaller, #11 or 
#2 | 

While there are always many organizational difficulties 
associated with writing and debugging a large program, the situation 
discussed above presented one of the most challenging theoretical 
problemsto resolve. The scheme finally chosen is stated below: 


If a combination's relative cost is greater than the current upper 
bound and the combination's identifying number is equal to 


i 

1) 2 , where i is any nonnegative integer, terminate genera- 
tion of further combinations (because all their relative 
costs will exceed the current upper bound), or 

ai ] : : 

2) 2° 23, where i and j are any nonnegative integers such 
that j < i, skip to combination number 21+! and continue 
generating, or 

3) otherwise, continue generating. 

In most cases, the above procedure greatly reduces the number of 
combinations necessary to generate. It could be expanded by saying 
- pal | Lue — 
that if the combination number is equal to 2 +2°+2 , where i, j, and 

; om ; i ae 
k are nonnegative integers such that k<j<i, skip to number 2 +2 ; 
and so forth. Except for very large problems, however, such added 
checks would probably consume more computer time than they would save. 

Another facet of the problem of trying to generalize the 

relationships among relative costs is presented when a combination is 
found constraint feasible. To illustrate the situation, suppose 
combination number 11 is found to provide a relative cost less than the 


current upper bound and is also found constraint feasible. If #11 is 


sent to SOLCHK at this point and found to provide a tour, there can be 
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no certainty that it is a least-cost tour among the current set of all 
possible combinations, This is because #12, as yet unchecked, may 
provide an even smaller relative cost, may be constraint feasible, and 
may also provide a tour. 

The technique devised to resolve the above difficulty was to 
store the identifying number and relative cost of each combination 
found constraint feasible. With each addition to storage, the numbers 
and relative costs are reranked in ascending order of relative costs. 
Finally, when all necessary combinations have been generated and 
checked, the constraint feasible solutions are sent to SOLCHK in 
ascending order of relative costs, thereby insuring the first tour 
found has a cost less than or equal to any tour the remaining solutions 
in storage might provide. As soon as a tour is found, control is 
returned to the first address indicated in the parenthesis of the 
calling statement in Main, If none of the solutions in storage provides 
a tour, control is returned to the statement in Main following the 
calling statement, 

Finally, a word about the regeneration of combinations required 
by the above procedure, The decimal identifying number in storage is 
converted to binary through the standard technique of successive 
divisions by 2. Remainders become the binary digits from right to left. 
The binary number, then, is the original representation of the combina- 
tion represented temporarily in storage by the decimal identifying 
number. Using the binary number, the combination of interest is re- 
constructed and sent to SOLCHK along with its associated relative 


cost. 


by 


C. LIMITATIONS OF THE PROGRAM 

The storage areas of the program are set up to deal with problems 
smaller than 41 cities, but these areas can be expanded easily if 
necessary. The four subroutines used to provide an optimal assignment 
solution are programmed to terminate execution in the unlikely event 
that more than 299 iterations are required. The only other significant 
limitations of the program are found in Subroutine LINCOM, but one of 
these presents a serious restriction. 

In Subroutine LINCOM, up to 1000 constraint feasible solutions can 
be stored at one time, but this does not appear to present an important 
limitation. In any case, storage space for CSTLC and NLC in Common 
Block TSBLK can be increased considerably if necessary. The most 
serious restriction and the one hardest to correct is the fact that 
Subroutine LINCOM can generate all the possible combinations for no 
more than 32 nonbasic variables. This means that if a tour is not 
found upon setting the 33 ranked, nonbasic variable equal to one, 
then execution on this problem must terminate. The reason for this 
Situation is that each combination must have a unique identifying number 
associated with it in LINCOM and the largest integer number that can 
be expressed on the IBM 360 is a There are 974.4 possible 
combinations of the S388 variable with the 32 preceding it, hence, 
the problem. As an example of the gravity of the restriction, note 
that a forty-city problem presents 1481 nonbasic variables for 
consideration. 

The perplexity above is not insurmountable, but it was not re- 
solved in this study. The possibility of representing the identifying 


number in floating-point mode is precluded because a form of the 
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number is used as a DO Loop parameter. (See the variable KT in LINCOM) 
It might be possible to represent the number as an element of a vector. 
For example, each of the first g31.y numbers might be prefixed with a 
1, the second set of et numbers with a 2, and so forth. The DO 

Loop using the number as a parameter could, in like manner, be 


expanded into a double DO Loop. 
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IV. COMPUTATIONAL EXPERIENCE 


Unfortunately, time constraints on this study did not permit 
extensive testing of the computer program in its final form. Enough 
was accomplished, however, to indicate several important conclusions. 

Appendix B presents a partial facsimile of the program output for 
the same ten-city problem that is discussed in References 8 and 9. 
Eastman's algorithm required the solution of eleven assignment prob- 
lems before the optimal tour was obtained. Little's method required 
31 branches to arrive at the desired solution. Appendix B indicates 
the programmed algorithm under study took only 1.50 seconds after the 
initial assignment solution to obtain the optimal tour. The first 
tour was found after trying combinations of the first 23 ranked, non- 
basic variables and this tour turned out to be the optimal one, 

City-to-city costs were generated by Subroutine RAND for 17 
twenty-city problems. Fourteen of these problems were solved by the 
program in a mean-time of OS seconds and after trying an average of 15 
nonbasic variables for each problem. Execution was terminated on the 
other three problems because a tour had not been found within the first 
33 nonbasic variables. 

Ten randomly generated thirty-city problems were attempted and 
nine were solved in an average time of 56 seconds. The mean number 
of nonbasic variables required for testing in these nine problems was 
MO. 

Finally, Subroutine RAND was used to create three forty-city 
problems. The first ran .07 seconds after the assignment solution to 


discover this was also an optimal tour. The second problem ran 549 
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seconds before terminating after trial of the 3374 nonbasic variable 
failed to provide a tour. The third problem was solved in 187 
seconds after trying only four nonbasic variables. 

In view of the computational experience cited above, one can 
conclude that the algorithm works and that the program works with one 
exception. Of the 31 problems tried, five could not be solved because 
tours had not been found upon setting the 33°4 nonbasic variables 
equal to one. If the program is to be of general use, this somewhat 
arbitrary restriction in Subroutine LINCOM must be removed. The 
limitation could be withdrawn if a way is found to obviate the 
identifying numbers presently associated with each possible combina- 
tion, Since oe aa is the largest integer that can be represented on 
the IBM System/360, none beyond the age ranked, nonbasic variable can 
be tried in combination with its preceding variables, 

Several steps can be taken to increase the execution speed of the 
program. All write statements, except those in Subroutine SOLCHK and 
the various termination messages throughout the program, can be replaced 
with CONTINUE statements. In subroutine REPART there is no particular 
reason why relative costs are recomputed and reranked after each 
artificial replacement. These two chores could be postponed until all 
artificials have been replaced. Since a replacement variable enters 
the basis at the zero level, it is not necessary that its relative 
cost be the least possible. 

In conclusion, a possible improvement to the algorithm, itself, is 
suggested. The algorithm is very efficient if a tour is found within 
the first few nonbasic variables tried, otherwise much time is consumed 


rechecking combinations already considered against a new upper bound, 
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and perhaps again against an even newer upper bound, and so forth. 
This phenomenon can be observed in the example output displayed in 
Appendix B. The algorithm might be improved as follows: 1) obtain 
a relatively low-cost tour by linking together, on some reasonable 
fashion, the subtours appearing in the optimal assignment solution, 
2) compute ® for this totir, 3) Ussing = as an upper bound, begin 
with the last, ranked, nonbasic variable whose relative cost is less 
than es and proceed backward just as in the original algorithm after 
the first tour is found, In this manner, it is not necessary to 
proceed forward through the ranked, nonbasic variables looking for 
that first tour, and it is necessary to work backwards through the 


variables only once. 
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Q STARTS CLOCK; 


SUBROUTINE TIMEIT(N, TIME) 
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APPENDIX B COMPUTER OUTPUT FOR A SELECTED PROBLEM 


2 
COST MATRIX 


1 2 3 4, 5 6 7 8 9 10 

l 299 24 18 tape 31 Ibe, 36 JAS 30 26 
2 1 999 1 Zi 26 SZ 25 J] 28 18 
o 22 ZS 999 Zo 16 Zo zy 18 16s, 
4 24 31 18 392 19 1S} 28 9 19 27 
5 18 18 34 20 229 ay 24 15 25 8 
6 24 1 iy 15 10 999 11 16 Zi 31 
y 28 es) 2] Bhs) |e 18 99% 21 21 19 
8 1g 24 18 163 ie ZZ Zs 399 mS, 24 
3 17 21 18 24 ZF 24 34 or 999 18 
10 18 19 20 16 Zo eye 18 31 23 999 


OPTIMAL ASSIGNMENT SOLUTION 


BASIS B 
82 0 
lik 
ise 
Shejmed: 
Cael 
29 

haem 
oa 
6a 
3) le, 
74 0 
62.sal 
SoPL 
0 

0 

0 

1 

0 

0 


Z 140 


Va : 
Alternate solution procedures for this ten-city problem are 
discussed in References 8 and 9. 
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TRY X( 8) = 1 REL COST = 0 

TRY X(97) = 1 REL COST = 1 

TRY THE (2** 1) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 2 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(68) = 1 REL COST = 2 

TRY THE (2** 2) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 2 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(99) = 1 REL COST = 2 

TRY THE (2** 3) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST ,LE. 2 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(20) = 1 REL COST = 2 

TRY THE (2** 4) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 2 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(84) = 1 REL COST = 2 

TRY THE (2** 5) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


oy 
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~ 


re 


j— 


Lo 


lo 


Lo 
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RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST ,LE. 3 = 
NR FOUND CONSTRAINT REASIBLE = 


TRY X(58) = 1 REL COST = 3 

TRY THE (2** 6) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X( 2) =m GREL COST =@ 

TRY THE (2** 7) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X( 9) = LoeREEe Gest es 

TRY THE (2** 8) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(19) = 1 REL COST = 3 

TRY THE (2%* 9) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(13) = 1 REL COST = 3 

TRY THE (2**10) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COSi=eyree3 | = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(52) = 1 REL COST = 3 

TRY THE (2**11) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 3 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(65) = 1 REL COST = 3 

TRY THE (2%*12) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(52) AND 

TRY THE (2**11) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


DZ 


OW & oo ee r te on aa a 


OW & 


RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(19) AND 

TRY THE (2%* 9) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 2) AND 

TRY THE (2%* 7) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 
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RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(28) = 1 REL®@OST = 4 

TRY THE (2%**13) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL GOSH .LEN® 49 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(96) = 1 REL COST = 4 

TRY THE (2%*14) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(17) = 1 REL COST = 4 

TRY THE (2**15) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 4 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(92) = 1 REL COST = 4 

TRY THE (2**16) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST ,LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(17) AND 

TRY THE (2**15) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(96) AND 

TRY THE (2**14) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(28) AND 

TRY THE (2**13) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST «LE 5 & 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(65) AND 

TRY THE (2**12) - 1 EN GoM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 
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OW 
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RETURN TO X(52) AND 

TRY THE (2%**11) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(19) AND 

TRY THE (2** 9) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 2) AND 

TRY THE (2** 7) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(20) AND 

TRY THE (2%* 4) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 
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RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(86) = 1 REL COST = 5 

TRY THE (2**17) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(54) = 1 REL COST = 5 

TRY THE (2**18) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(66) = 1 REL COST = 5 
TRY THE (2**19) - 1 LN COM 
_ NR LN COM NECESSARY TO GENERATE 
. NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(71) = 1 REL COST = 5 

TRY THE (2**20) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(25) = 1 REL COST = 5 

TRY THE (2**21) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 5 = 
NR FOUND CONSTRAINT FEASIBLE = 


TRY X(85) = 1 REL COST = 5 
TRY THE (2**22), - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(25) AND 

TRY THE (2**21) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 
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RETURN TO X(71) AND 

TRY THE (2**20) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(66) AND 

TRY THE (2**19) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(54) AND 

TRY THE (2**18) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(86) AND 

TRY THE (2**17) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(92) AND 

TRY THE (2**16) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(17) AND 

TRY THE (2**15) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(96) AND 

TRY THE (2**14) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


A TRAVELING SALESMAN SOLUTION IS -- 


X( 2, 1) = 
X( 8, 5) = 
X( 4, 8) = 
X( 5,10) = 
FD i 
X( 1, 3) = 
X( 7, 2) = 
Ge a 
X(10, 6) = 
1 i) 


Ons) RO 


ws we Le 


SS ee ee eR ee eS 


Z = 146 
oy 
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me WwW 


16 
ll 


16 
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RETURN TO X(28) AND 

TRY THE (2**13) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(65) AND 

TRY THE (2**12) - 1 LN COM 
NR LN GOM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(52) AND 

TRY THE (2**11), - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL GOST .LEY” 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(19) AND 

TRY THE (2** 9) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X( 2) AND 

TRY THE (2** 7) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND ‘CONSTRAINT FEASIBLE = 


RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN GOM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE, 6 = 
NR FOUND CONSTRAINT FEASIBLE = 
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RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST .LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 


RETURN TO X(97) AND 
MRY THE {2%% 1) -9 1) EN Com 
NR LN COM NECESSARY TO GENERATE 
NR WITH REL COST. LE. 6 = 
NR FOUND CONSTRAINT FEASIBLE = 
*kZ = 146 IS OPTIMAL ** 


ALGORITHM TIME = 1.5043 SECONDS 


oy) 


US, 
13 


“JI ™“ 


O W 


te 


LIST OF REFERENCES 
Dantzig, George B., Linear Programming and Extensions, p. 545-547, 
2d ed., Princeton University Press, 1966, 
5 woneeOS:. 
» Chapters 14 and 15. 


Ford, L. R., Jr., and Fulkerson, D. R., Flows in Networks, 3d ed., 
Chapter III, Princeton University Press, 1967. 


Hadley, G., Linear Programming, 2d ed., Chapter 9, Addison-Wesley 
Publishing Company, Inc., 1962. 


Hadley, G., Nonlinear and Dynamic Programming, p. 267-269, Addison- 
Wesley Publishing Company, Inc., 1964. 


IBM Systems Reference Library, IBM System/360, Fortran IV Language, 
File No. S360-25, Form C28-6515-5, 1965, 1966. 


Lawler, E. L. and Wood, D. E., "Branch-and-Bound Methods: A 
Survey,'' Operations Research, v. 14, p. 707-709, July-August, 
1966. 


, p. 709-711. 


60 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Civil Schools Branch 
Department of the Army 
Washington, D. C, 20315 


Professor Harold Greenberg, Code 55Gd 
Department of Operations Analysis 
Naval Postgraduate School 

Monterey, California 93940 


Captain William W. White, USA 


2040 Hanover Circle 
Beaumont, Texas 77706 


61 


No. Copies 


20 





Security Classification 









DOCUMENT CONTROL DATA-R&D 


(Security classification of title, body of abstract and indexing annotation must be entered when the overal! report Is classified 


24. REPORT SECURITY CLASSIFICATION 
Unclassified 
eee: 






1. ORIGINATING ACTIVITY (Corporate author) 


Naval Postgraduate School 
Monterey, California 93940 









3. REPORT TITLE 






A Traveling Salesman Algorithm 






4. OESCRIPTIVE NOTES (Type of report and, inclusive dates) 


Master's Thesis; October 1969 


§. AUTHOR(S) (First name, middle initial, last name) 


William Willerson White 


6. REPORT OATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
October 1969 62 9 





6a. CONTRACT OR GRANT NO. 9a. ORIGINATOR'S REPORT NUMBER(S) 





6b. PROJECT NO. 


95. OTHER REPORT NOI(S) (Any other numbere thet may be aesigned 
thia report) 


10. DISTRIBUTION STATEMENT 


This document has been approved for public release and sale; its distribution 
is unlimited. 
tt. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


13. ABSTRACT 


A possible approach to solving the traveling-salesman problem was suggested by 
Professor Harold Greenberg. This study developed an algorithm based on his suggestion 
and then programmed the algorithm for the IBM System/360. While the program solves 
problems of forty cities or less, it has a significant limitation, Execution is 
terminated on a problem if a solution is not found early enough in the trial-and- 
error process of the algorithm. The solution procedure developed formulates the 
salesman's problem as an assignment problem, obtains an optimal assignment solution, 
and then manipulates vectors in the final simplex tableau until an assignment 
solution is obtained that also satisfies the additional traveling-salesman 
constraints. Background to the problem is given, the algorithm is developed and 
stated, the computer program is described and critiqued, highlights of computational 
experience with the program are presented, and, finally, some conclusions and 
recommendations are made, 


DD ar.1473 Par . 
S/N 0101-807-6811 Security Classification | | 


4-31408 


~ Security Classification 


KEY WOROS 


Algorithm 


Linear Programming 


DD ™.1473 orc — 


S/N 0101-807-6823 


Security Classification A-31409 

















thesW55542 
A traveling salesman algorith 


726 
DUDLEY KNOX LIBRARY 





